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Abstract 

Previous papers dealt with the quantization of the Lemaitre-Tol- 
man-Bondi (LTB) model for vanishing cosmological constant A. Here 
we extend the analysis to the case A > 0. Our main goal is to present 
solutions of the Wheeler-DeWitt equation, to give their interpretation, 
and to derive Hawking radiation from them. We restrict ourselves to 
a discussion of those points that are different from the A = 0-case. 
These have mainly to do with the occurrence of two horizons. 
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1 Introduction 



Part of the research on quantum gravity consists in the attempt to get a 
possible insight into the final theory from concrete models. The simplest 
models deal with spatially homogeneous cosmological metrics and are typ- 
ically used in quantum cosmology. Less simple are models with spherical 
symmetry. They are used in cosmology, too, but mainly in the description 
of black holes. 

Here we deal with a particular spherically-symmetric model, the Le- 
maitre-Tolman-Bondi (LTB) model. It describes self-gravitating inhomo- 
geneous dust and is well understood at the classical level where it is mostly 
used in cosmology [1] . The quantization of this model was attempted in both 
quantum geometrodynamics and loop quantum gravity. Here we restrict 
ourselves to quantum geometrodynamics, in which the central equations are 
the Wheeler-DeWitt equation and the momentum (diffeomorphism) con- 
straints [Ills]. 

Following the application of canonical geometrodynamics to the Schwarz- 
schild black hole [5] and to Reissner-Nordstrom-anti-de Sitter black holes 
[S], the canonical formalism for the LTB model was developed in ^ and then 
applied to quantization in a series of papers, see El |9] and the references 
therein. While it was not possible to construct the quantum theory for 
the LTB model in a rigorous way, partial progress was achieved. Among 
the results was the recovery of Hawking radiation plus greybody corrections 
from exact solutions to the Wheeler-DeWitt equation and the momentum 
constraints \E, 15]. They were found for the case of vanishing cosmological 
constant A where only one horizon (the black- hole horizon) is present. Here 
we extend the analysis to the case A > 0, where two horizons (the black- 
hole and the cosmological horizon) exist. We again manage to find exact 
quantum states from which the two different Hawking temperatures from 
the two horizons can be recovered. 

Although we are mainly interested in understanding quantum gravity, 
we want to emphasize that the case A > fits very well current cosmological 
observations [10]. The analysis for A < 0, where only one horizon is present, 
was performed in [llj . 

In our present paper we shall focus on the main differences to the case 
of vanishing A. For more details on the general formalism we refer to the 
earlier papers cited above. The LTB model was also addressed from the 
perspective of loop quantum gravity [T2] . 
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2 The classical LTB model with positive cosmo- 
logical constant 



Since we deal with dust, the energy-momentum tensor is given by T^^, = 
e(r, p)u^Uu-, where = ■u'^(r, p) is the four-velocity vector of a dust particle 
with proper time r; the parameter p labels the various shells which together 
form the dust cloud. The line element for the LTB spacetime in comoving 
synchronous coordinates is given by 



d.^ = -dr^ + T^^^dp2 + R\p, r){de' + sin^ 9d<p') . (1) 
1 + 2E[p) 

Inserting this expression into the Einstein field equations leads to 



and 



2 F Ai?2 



{drRY = - + ^ + 2E = l-T+2E , (3) 

where F{p) = 2GM{p) is a non-negative function with the dimension of a 
length, and 

The case of collapse is described by drR{j, p) < 0. We set c = 1 throughout. 

There still exists the freedom to rescale the shell index p. This freedom 
can be removed by demanding 

R{0,p)=p, (5) 

so that for r = the label coordinate p is equal to the curvature radius 
R. Now we can express the functions F{p) and E{p) in terms of the energy 
density e at r = 0. From ^ and ([3]) one gets 

F{p) = 87tG r e{0,p)p'dp , (6) 
Jo 

E{p) = l[drRir = 0,p)f-^-^. (7) 

The interpretation of these quantities is that F{p)/2G = M{p) is the active 
gravitating mass within a p = constant shell, while E{p) is the total energy 
per unit mass within the same shell. The marginally bound models are 
defined by E{p) = 0. Here we consider the general case, which includes the 
non-marginal case defined by E(p) ^ 0. 



2 



In order to derive Hawking radiation, we shall consider the following sit- 
uation. We assume the presence of a black hole with mass Mq surrounded 
by a gravitating dust cloud. We assume, moreover, that the total energy of 
the cloud is small compared to the mass of the black hole. After quantiza- 
tion, we shall find exact quantum states for the dust cloud from which the 
Hawking temperatures will be found. The dust cloud will play the role of 
the quantum fields usually employed in the derivation of Hawking radiation. 

If the dust cloud is of finite extension, the metric outside of it is of the 
Schwarzschild-de Sitter (SdS) or Kottler formH 

dT^ 

-1 

di?2 + R^dn^, (8) 

where M is the mass of the black hole plus the mass contribution of the 
cloud. For small R, SdS space approximates Schwarzschild space, while 
for large R and positive A it approximates de Sitter space. The black-hole 
horizon i?h and the cosmological horizon Rc are two of the three zeros of 
qtt in &, see, for example, [H] and [TSj. They are explicitly given by 

where ^"^ = 3GM^/A and ^ = cos cos ^]). The third zero value, Rn, 
is negative, R^ = —R^ — Rc = —QGMl^ and possesses no obvious meaningo 
The LTB solution ([1]) must then be matched at the dust boundary to the 
SdS solution. 

In the canonical formalism, the general ansatz for a spherically-symmetric 
line element reads 

ds^ = -iV^dt^ + l2 _ ^r^^^2 ^ ^2^^2 ^ ^^q^ 

where is the lapse function, A^** is the only component of the shift vector 
that survives the symmetry reduction, and dJl^ = d^^ -|- sin^ Odcf)^ is the line 
element on the unit two-sphere. 

^See, for example, [T^ for a detailed discussion of this metric. 

*It must be assumed here that the black-hole horizon is smaller than the cosmological 
horizon, i?h < Rc, which means that < 1. This is well satisfied in our Universe: 
Inserting for M the mass of the supermassive black hole in the quasar OJ287, which is 
presently the biggest known supermassive black hole, with M ~ 1.8 x IO^^Mq [16], and 
for A the value following from ~ 0.726 [ID], we get ^ 9 x 10~". The value 1 = 1 
would correspond to the Nariai metric for which both areas are the same but are separated 
by a finite proper distance. (See [T3 for a discussion of the Nariai metric and a reprint of 
Nariai's original papers.) 



ds 



2 



2GM AR^ 



+ 1 



R 3 
2GM Ai?2 
~R 3~ 



Rh = 
Rc = 



3 



Inserting the ansatz ()10p into the ADM form of the Einstein-Hilbert 
action, we obtain the gravitational part of the action, 



oo 



S^= dt dr[PLL + PRR-NrL^-N''nf]+Sdj^, (11 





where the Hamiltonian and the diffeomorphism (momentum) constraint are 
given by 

^ --^i-?r-^j + Gr2 - 2L + V W ' ^^^^ 

nf = r'Pr - LP'l , (13) 

respectively (a prime denotes a derivative with respect to r). The boundary 
action Sq^ will be discussed below. 

The total action is the sum of pi|) and an action S'^ describing the 
dust. The canonical formalism for the latter was developed in [18] . The 
dust action reads 

3"^ = j dt dr [Prf - NW^ - N^'T-if^ , (14) 
where the contributions to the Hamiltonian and momentum constraints are 

n'^ = P^Jl + — and nf = T'Pr. (15) 

It is due to this relatively simple form of the dust action that exact quantum 
states can be found below. The full constraints are then the sum of the 
gravitational and the matter (dust) constraints. 

The general ansatz of the metric (fTOj) should, of course, correspond to 
the LTB metric ([1]). This leads to a couple of relations that are analogous 
to those derived in j7|, except that is now given by One can express, 
in particular, the variables E, F, and r, which completely characterize the 
collapse of the dust cloud, in terms of the canonical variables. 

What about the boundary action Sgj] in ()lip ? In principle, we have to 
consider two boundaries. The first boundary is at the centre of symmetry 
for the dust cloud, the second either at spatial infinity or at another chosen 
end of the spatial hypersurfaces. As for the inner boundary, the effect of 
a positive cosmological constant compared to a vanishing or negative A is 
there negligible, because its effects are felt only at cosmological distances. 
We can thus employ the boundary action used previously, cf. [9j. It was 
found that the following term must be added to compensate a corresponding 
term that arises from variation of the action and partial integration: 

5as = [ dt No{t)Moit) , (16) 
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where Nq is the lapse function at r = 0, and Mq = F{0)/2G is the mass of 
the black hole at the centre. If this compensation were not made, one would 
conclude that a variation of A^o leads to vanishing black-hole mass, which 
would certainly not be desirable. Such arguments were already used in the 
geometrodynamics of the pure Schwarzschild black hole [1] . 

The choice of the second boundary is less obvious. Since the spacetime 
here is asymptotically de Sitter, the spatial sections are no longer asymp- 
totically flat. This leads basically to two options. One option is to choose a 
spacelike hypersurface that approaches the spacelike infinity in SdS space- 
time. Here one must be careful and avoid r to become timelike. This is, for 
example, achieved if we employ coordinates such as the Painleve-Gullstrand 
coordinates used in [19] which are, in fact, close to the coordinates that we 
use below to simplify the constraints. Another option is to use the cosmolog- 
ical horizon as the second boundary. This would correspond to the choice 
frequently made in black-hole papers where the bifurcation sphere of the 
horizon is employed as the inner boundary, see, for example, |5j and [20| . It 
will also be our choice here. 

We thus consider the case where the cosmological horizon is the second 
boundary. The situation is then analogous to the cases considered in [5] and 
|2Uj . We choose r = to represent the centre of the cloud and let r ^ oo 
at the cosmological horizon. This is possible if we smoothly match r to 
the tortoise coordinate i?*, defined by dR^: = J^^^dr, at the cosmological 
horizon. The time parameter t is chosen equal to the Killing time T. Let us 
now turn to the fall-off conditions. As for L, we choose 

i(,0 = i£«+o(^) . (17) 

L tends to zero near the cosmological horizon, where r goes to infinity. 
This is required, since it is equal to the tortoise coordinate there. It is 
proportional to which guarantees that the distance to the horizon will 
be finite. As for R, we choose 

= + ^+0(l) . (18) 

For the remaining variables we impose 

-('•r^-^^O^r^, N'.0(^). (19) 

and 

These fall-offs are consistent with the equations of motion and the constraint. 
They also guarantee that the boundary term at the cosmological horizon 
vanishes. 
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One could now start with the quantization of ()12p . (|13p . and (|14p . How- 
ever, this turns out to be too comphcated. One thus performs first a classical 
simplification of the constraints in order to render them manageable. This 
is achieved by introducing new variables and appropriate canonical trans- 
formations [U [7] . In this way, the variable T = F' and its momentum Pp 
are introduced as new canonical variables; this choice is also convenient 
because one thereby absorbs a certain boundary term [7]. Using then the 
momentum constraints to eliminate Pr and squaring the ensuing Hamilton 
constraint, one arrives at the following final form of the constraints: The 
full Hamiltonian constraint reads 

W = P2+^P2__^^0^ (21) 

while the full momentum constraint reads 

Ur = R'Pr - rPr + r'Pr w . (22) 

This is the form of the constraints suitable for quantization. Strictly speak- 
ing, re-defining the constraints is part of the very definition of quantization. 
There are many cases where the classical constraints are transformed into a 
manageable form by canonical transformations, see, for example, |21j . 

We note that because of the involved squaring the new Hamiltonian con- 
straint (j2ip has acquired the dimension (mass/length)^. We emphasize that 
the only differences compared to the earlier papers dealing with non-positive 
A lie so far in the definition of J- according to We also emphasize that 
through these manipulations the kinetic term in the Hamiltonian constraint 
is no longer hyperbolic, in contrast to the original form (jl2p . More precisely, 
it is hyperbolic inside the black-hole horizon and outside the cosmological 
horizon, and it is elliptic between the horizonsH A similar observation was 
made for the Reissner-Nordstrom black hole already in j22]FI This gives 
further support for choosing a hypersurface that extends from one horizon 
to the other. We finally note that the gravitational constant G occurs in 
()2ip in the same way as in the usual Hamiltonian constraint without mat- 
ter; this may be of relevance for semiclassical approximation schemes in the 
quantum theory 0. 

The SdS spacetime is static only in between the two horizons. It thus 
makes sense to talk about a Killing time, T, only in this region. In previous 

^One might think that this change of sign is related to the fact that the metric in the 
interior of the horizon can become "cosmological", as it happens e.g. for the Schwarzschild 
black hole where the inside metric is of the Kantowski-Sachs form. This is, however, not 
the case. The equation .F = is a condition for an apparent horizon, and the sign change 
happens also for the non-static case where the interior is not "cosmological" . 

^It is interesting to note that the Einstein equations for regular axisymmetric and 
stationary black holes surrounded by matter are also elliptic in the exterior and hyperbolic 
in the interior of the hole |23) . 
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papers it was shown that T = 2Pr and that 



T = 



ar ± 



dR 



(23) 



where a = l/Vl + 2£'; the plus sign in ()23p holds for an expanding, the 
minus sign for a collapsing dust cloud. (This interpretation comes from the 
similarity of the coordinates (r, R) with the Painleve-Gullstrand coordinates 
for a Schwarzschild black hole [21] •) It is straightforward to show that the 
same relation also holds in our case, taking into account the new definition 
for Strictly speaking, one can talk about a Killing time only in the 
absence of the dust cloud, where the spacetime is of SdS form. However, as 
in the previous papers, one can approximately continue to interpret T as a 
Killing time, because we consider the dust only as a small perturbation to 
the SdS spacetime. 

3 Quantization 

We now turn to quantization. Applying the formal Dirac procedure, we 
replace the momenta by functional derivatives with respect to the corre- 
sponding configuration variables. From the Hamiltonian constraint ()2ip we 
get the Wheeler-DeWitt equation. 



Here, A and B are smooth functions of R and F, which encapsulate the 
factor-ordering ambiguities. The factor-ordering problem is a fundamental 
problem that can only be dealt with after a suitable regularization procedure 
has been invoked. A general treatment is beyond the scope of our paper. 
To implement the factor-ordering ambiguities at a formal level, we have 
introduced formal factors of (5(0) into the equation. The reason is that the 
equation can then be put onto a lattice and the continuum limit can be 
performed [7]. This corresponds to a particular regularization procedure 
that allows the presentation of the solutions below. 

From (|22p one gets the quantum momentum constraint 



Both equations, (p^ and (f25|) . can be put on the lattice and a special class 
of exact solutions can be found In the continuum limit, this class of 




(24) 




(25) 
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solutions reads 



^[T{r),R{r),T{r)] 
i 



exp lb 



2Gh 



dr r 



ar + / di? 



(26) 



(In the earlier papers, a more general solution with an unspecified constant 
b was used [7j; for simplicity, we set here 6 = 0.) The momentum constraint 
([25|) is solved by any state of the form exp (i/{2GH) f dr TW), where W is 
a smooth function of R and r. The Wheeler-DeWitt equation, however, 
is only solved for a particular factor ordering. Generalizing the treatment 
presented in [7] to A 7^ 0, one must choose in B = and 



A 



F 



2AR\ 2 - a'^T 



(27) 



This means that A does not only enter the definition of T ^ but also the 
chosen factor ordering. Exact solutions for the quantum states can only be 
found for this special A-dependent choice of A. This is certainly a weak 
point of our approach, because there should be a fixed factor ordering at the 
fundamental level. Nevertheless, this choice for A gives us the opportunity to 
have exact solutions at our disposal. Moreover, at the semiclassical level of 
WKB states, the factor ordering is irrelevant, and the calculation of Hawking 
radiation below is an effect that occurs at the WKB level |25J. In fact, one 
can use a semiclassical expansion scheme instead of the exact solutions to 
derive Hawking radiation and quantum gravitational corrections [26j. 



4 Hawking radiation 



Let us now turn to the calculation of Hawking radiation. The idea is to 
remain as close as possible to Hawking's original derivation in [27], but to 
transfer this idea to solutions of the Wheeler-DeWitt equation and the mo- 
mentum constraints, which are full quantum gravitational states. Because 
the Wheeler-DeWitt equation (j24|) contains second derivatives with respect 
to the dust time r (which is a consequence of squaring the Hamiltonian at 
the classical level), we can define states of positive and negative frequency 
with respect to r. This is crucial for our discussion. The situation is thus 
analogous to the Klein-Gordon equation, where positive frequency is defined 
by a minus sign in the phase, exp(— ii?t/^), and correspondingly negative 
frequency by a positive sign. 

For the state with positive frequency we thus get from ([26]) . 



[r(r),i?(r),r(r)] 
/ i 



exp 



dr r 



ar + I dR — 



(28) 
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while for the state with negative energy we have 



^- [T{r),R{r),T{r)] 



ar + I dR — 



(29) 



In Hawking's calculation in |27j . the Bogolyubov coefficient /5, which gives 
the negative-frequency part of the time-developed original state with posi- 
tive frequency, plays a central role. Here, the role of the quantum matter 
field in [27] is played by the dust. We thus define Hawking radiation for 
the black-hole the overlap between an "outgoing dust state" with 

negative energy and an "ingoing dust state" with positive energy. Since the 
interpretation of these states is made with respect to an observer in the SdS 
spacetime using the Killing time T, we have to substitute the dust time by 
T according to ()23p . Taking the corresponding signs in (|23p into account, 
we have then to calculate the overlap between 



*e" = exp ( ^ / dr rr ) , (30) 



which is actually independent of R, and 



where the index e(c) denotes "expanding" ("collapsing"). 

For the cosmological horizon, the situation is just the opposite because 
the Hawking radiation is now incoming from the horizon. We thus have to 
calculate the overlap between the "ingoing negative-energy state" 



1 



R 



and the "outgoing positive-energy state" 



In order to calculate the overlap between the quantum states ()30p and 
([3T]) . or between ([32]) and ([33]) . we need a Hilbert space. As far as full 
quantum gravity is concerned, this is an open issue [2]. In the present 
case, however, a natural candidate is present [251 E]- We can choose the 
standard Schrodinger inner product with respect to constant time and R 
as the integration variable. The integration range is supposed to run from 
the black-hole horizon, i?h, to the cosmological horizon, R^. This is further 
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motivated by the fact that our hypersurface hes between the horizons; as 
we have already noted, in this region the kinetic term in (j2ip is ehiptic. 

Addressing first the cosmological horizon, we have for the desired over- 
lap between ingoing negative-energy and outgoing positive-energy states the 
expression 



Here, guji is the i2i?-component of the DeWitt metric, as it can be read 
off (j2ip where its inverse is the prefactor of the term P^, and we thus have 



There is, however, one further point to consider. Since we are interested 
in an observer far away from the dust cloud, for whom the Killing time 
T according to ([23]) is the appropriate time coordinate, we have already 
rewritten the states in terms of T. Consequently, as in [251 E]) we have to 
evaluate also the component g^R in the new coordinate system (T, R) instead 
of {t,R). Using ([25]) . one gets ^qrr = {aT)~^ . This is the expression to 
be used in the calculation of the Bogolyubov coefficient. 

For the concrete calculation we shall write the full states as a product 
of single-shell states where the radial variable r is assumed to consist of 
discrete points separated by a distance a. (The continuum limit is obtained 
for a — > 0.) As in [8], the Bogolyubov coefficient /3 is calculated for each 
shell separately. In the discrete case, we replace T by the dimensionless 
variable 2uj and indicate the dependence on tu by an index. (The factor 2 
is motivated by the fact that T = 2M'.) We omit the shell index and write 
the corresponding wave functions as -i/'i^ (T, i?) . From (j34p we then define /3 
to read 



In the earlier papers [8] and [9j we have chosen a particular normalization for 
/?, which was motivated by the normalization of states for the Klein-Gordon 
equation. Since the full normalization of solutions to the Wheeler-DeWitt 
equation is anyway not known, we shall leave the question of normalization 
open here and define /? directly by ([35]) . 

Employing now the one-shell contributions of (|32p and (|33p and inserting 

''Note that tliis measure is not the measure that would render the Wheeler-DeWitt 
equation (I24|l with the factor ordering ()27|1 Hermitean [9]. Our choice of the measure in 
(|34[) would make the _R-part of (|24[) Hermitean if the Laplace-Beltrami factor ordering 
were used. 




(34) 



9rr = :f-'E 




(35) 
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them into ()35p . one gets the following expression: 



ifjT, 



Vl + 2E exp ( - — [lj + w'] I X 



dR T ^ exp 



2icra; 



\.^]. (30) 



We shall first calculate the contribution to from the cosmological hori- 
zon. This will be achieved through an appropriate near-horizon approxima- 
tion, making the standard assumption that the integral is dominated by the 
close vicinity of the horizon. We start by writing in terms of partial 
fractions: 

1 1 



1 _ Z a_r2 



R 



R 3 

ABC 
+ — 7^ + 



R — Rc R — Rh R — 



(37) 



Now we introduce the variable shy Is = Rc — Rso that near the cosmological 
horizon |s| <C 1. Here, / = 1/\/A- 
First we consider the term 



-lds{Rc - Is) 



A 

-T- + 



B 



+ 



C 



\ 



1 



Rn — Is 



Is Rc — Rh — Is Rc — Rn — Is 

] ■ (38) 



1 



+ 



B 



+ 



C 



Is ^ Rc-Rh-ls ^ Rc-Rn-ls / 

After some calculation, neglecting terms with a positive power in s, we get: 



dsA/const - + Xa'^R'^ 



(39) 



with ^ = 2A'^lRc + 
yields 



+ ^r3! - a^MRc. Expanding ([MD for |s| <Cl 



dsW const + -^A^Rl 



is'- 

s 



AR, 



is 



2AR, 



+ 0{s' 



PS ARcln\s\ 

Inserting this in our integral (j36p and using 

AR ARc 



is 



2ARr 



(40) 



R — Rc 



Is 
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we obtain 



P^^, = -VTT2Eexp(^-i^{uj + uj')^AR, 
Jo 

To evaluate this integral, we use the formula |28j 



i,^ aujs 
ARr Gh 



(41) 



which is, in particular, applicable to the case p = and < Rez^ < 1. 
(We insert a small positive value for Rei^, which we let go to zero after the 
integration.) Then, 

(3^^i = -Vl + 2Eexp + ARc x 



Using 



r(-iy)r(iy) 



y sinh Try 
(with real y), we get 

nGhAR,il + 2E) 1 

for the absolute square of /?. (In spite of the minus sign, this is positive 
because A < 0, see below.) Next we want to calculate the particle creation 
rate, which means we have to evaluate the expression Sfc/3jfc/9*j,. Here this 
corresponds to the integral /q°° da;'|/3^^'p. This integral diverges and hence 
needs to be suitably regulated. We introduce a decay factor e~^' and carry 
out the integration. For the "in" particle number operator we then obtain 

(m iVin m) = — TTT^ • (44) 

Replacing auj by GAe, where Ae is the energy of a shell, we arrive at the 
final result 

m iVin m = ^ — . A ,^ , 45 

from which we read off 
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How can one interpret the prefactor in ()45p ? It should give the greybody 
factor, but its interpretation is comphcated by the fact that the fuh nor- 
mahzation of our wave functionals is open. We thus do not address it any 
further. 

Using ()37p . one can easily determine 
A 



A{Rc - Ri,){Rc - Rn) 
and thus finds 

= mR.-R.)iR.-R.) 

UttRc ^ ' 

This result coincides with Equation (2.15b) in [14]. 

Interestingly, our calculation yields a Planck spectrum with a tempera- 
ture that would be measured by a fictitious observer with qtt = 1- While 
such an observer could be easily realized in the Schwarzschild case (he would 
be situated at spatial infinity), this is not possible here. For the SdS case, 
the observer should be located between the horizons. Following [29j, we shall 
consider the preferred observer who is stationary at a position R^ where the 
black-hole attraction and the cosmological expansion cancel each other; such 
an observer is thus unaccelerated and moves on a geodesic. We hav^ 

which leads to the modified Hawking temperature 

= 1 M(fl,-fl,)(fi,-fl„) 

Vl - (9G2M2A)i/3 l2irB^ 

This also agrees with the corresponding expression for the temperature given 
in Equation (17) in [15] where the Hawking temperature is calculated from 
a tunnelling picture. 

A completely analogous calculation yields the Hawking temperature for 
the black-hole horizon. In fact, the corresponding overlap between ^I'" and 

leads to the same expression ()35p for (3. In contrast to above, however, 
the near-horizon approximation is now performed for the black-hole horizon. 
Now, the coefficient B in ()37p enters. It is given by 

B 



A(iih — Rc){Rh — Rn) 



*For the above example of the supermassive black hole in the quasar OJ287, one gets 
Ro ~ 260 kpc, which corresponds to roughly five times our distance to the Large Magel- 
lanic Cloud! 
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which follows from the expression for A by interchanging Rc and i2h- One 
thereby arrives at the expression 



hAiR,-R^){R^-R^) 

^^^^ = ' 

which is equal to Eq. (2.15a) in [H]. This expression together with (j47p 
are the main results of our paper, because they have been derived from 
candidates for exact quantum gravitational states. The temperature ()50p 
will be modified by the same factor as in ()49p if an observer is considered 
who follows the geodesic motion at Rq. 

The two temperatures (|47p and (|5Up arise from the same expression (|35p 
by two different near-horizon approximations. An exact evaluation of ()35p 
would yield the sum of both temperatures plus interference terms between 
the two types of Hawking radiation. In the extremal case of the Nariai 
metric, both temperatures are zero. Nevertheless, an exact evaluation of 
()35p could yield a small non-zero contribution which would be a genuine 
quantum gravitational effect. 

To summarize, we have shown that solutions to the Wheeler-DeWitt 
equation and the diffeomorphism constraints contain information about the 
two Hawking temperatures from the black-hole and the cosmological hori- 
zon. An interesting open question is the calculation of the entropy for the 
Schwarzschild-de Sitter case through a counting of microstates. Such a 
derivation was presented for the BTZ black hole [30] and the AdS black hole 
|31j . Whether it is also possible here is left for future publications. 
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